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The correspondence principle asserts that quantum mechanics resembles classical 
mechanics in the high-quantum- number limit. In the past few years many papers 
have been published on the extension of both quantum mechanics and classical me- 
chanics into the complex domain. However, the question of whether complex quan- 
tum mechanics resembles complex classical mechanics at high energy has not yet 
been studied. This paper introduces the concept of a local quantum probability den- 
sity p{z) in the complex plane. It is shown that there exist infinitely many complex 
contours C of infinite length on which p{z) dz is real and positive. Furthermore, the 
probability integral J c p{z) dz is finite. Demonstrating the existence of such contours 
is the essential element in establishing the correspondence between complex quan- 
tum and classical mechanics. The mathematics needed to analyze these contours is 
subtle and involves the use of asymptotics beyond all orders. 
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I. INTRODUCTION 



In conventional quantum mechanics, operators such as the Hamiltonian H and the posi- 
tion x are ordinarily taken to be real in the sense that they are Hermitian (H = W). The 
condition of Hermiticity allows the matrix elements of these operators to be complex, but 
guarantees that their eigenvalues are real. Similarly, in the study of classical mechanics the 
trajectories of particles are assumed to be real functions of time. However, in recent years 
both quantum mechanics and classical mechanics have been extended and generalized to the 
complex domain. In classical dynamical systems the complex as well as the real solutions to 
Hamilton's differential equations of motion have been studied [TUTU]. In this generalization 
of conventional classical mechanics, classical particles are not constrained to move along the 
real axis and may travel through the complex plane. In quantum mechanics the class of phys- 
ically allowed Hamiltonians has been broadened to include non-Hermitian PT-symmetric 
Hamiltonians in addition to Hermitian Hamiltonians, and wave functions (solutions to the 
Schrodinger equation) are treated as functions of complex coordinates [U I20H28] . Experi- 
mental observations of physical systems described by complex PT-symmetric Hamiltonians 
are now being reported [291436] . 
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FIG. 1: Graphical illustration of the correspondence principle. This figure compares the probability 
densities for a quantum and for a classical particle in a parabolic potential. The normalized 
quantum- mechanical probability density p(x) = \ipi(x)\ 2 for a quantum particle of energy = 9 
is plotted as a function of x (solid curve); p{x) is wavelike and exhibits four nodes. The normalized 
classical probability density for a particle of the same energy is also plotted (dotted curve); the 
classical probability density is not wavelike. Both probability densities are greatest near the turning 
points, which lie at ±3, because the speed of the particle is smallest in the vicinity of the turning 
points. Thus, the particle spends most of its time and is most likely to be found near the turning 
points. The classical probability density, which is the inverse of the speed s of the particle, is infinite 
at the turning points because s vanishes at the turning points, but this singularity is integrable. 

The relationship between quantum mechanics and classical mechanics is subtle. Quantum 
mechanics is essentially wavelike; probability amplitudes are described by a wave equation 
and physical observations involve such wavelike phenomena as interference patterns and 
nodes. In contrast, classical mechanics describes the motion of particles and exhibits none 
of these wavelike features. Nevertheless, there is a connection between quantum mechanics 
and classical mechanics, and according to Bohr's famous correspondence principle this subtle 
connection becomes more pronounced at high energy. 

Pauling and Wilson [37J give a simple pictorial explanation of the correspondence principle 
and this explanation is commonly repeated in other standard texts on quantum mechanics 
(see, for example, Ref. |38j). One compares the probability densities for a quantum and for 
a classical particle in a potential well. For the harmonic-oscillator Hamiltonian 

H = p 2 + x 2 (1) 

the quantum energy levels are E n = 2n + 1 {n = 0, 1, 2, 3, . . .). The fourth eigenfunction 
ip4(x) has four nodes, and the associated probability density p(x) = \ip4(x)\ 2 also has four 
nodes (see Fig. [I]). When n = 16, the quantum probability density has 16 nodes (see Fig. [2]). 
The probability density for a classical particle is inversely proportional to the speed of the 
particle and is thus a smooth nonoscillatory curve, while the quantum-mechanical probabil- 
ity density is oscillatory and has nodes. Figures [T] and [2] show that the classical probability at 
high energy is a local average of the quantum probability. (The nature of this averaging pro- 
cess can be studied by performing a semiclassical approximation of the quantum-mechanical 
wave function. However, we do not discuss semiclassical approximations here.) 

The objective of this paper is to extend and expand the work in a recent letter that exam- 
ines the relationship between complex classical mechanics and complex quantum mechanics 
[39] . In particular, we show how to generalize the elementary Pauling-and- Wilson picture of 



FIG. 2: Comparison of the normalized probability densities for a quantum particle and for a 
classical particle of energy E\& = 33 in a parabolic potential. The quantum probability density 
(solid curve) has 16 nodes and is wavelike while the classical probability density (dotted curve) is 
not wavelike; inside the potential well the quantum density oscillates about the classical density. 

the correspondence principle into the complex domain. We will see that in the complex do- 
main some of the distinctions between classical mechanics and quantum mechanics become 
less pronounced. For example, in complex classical mechanics, particle trajectories can enter 
classically forbidden regions and consequently exhibit tunneling-like phenomena [iOl W\\ . 

Figures [T] and [2] highlight an important difference between real classical systems and real 
quantum systems in the classically forbidden region: Outside the classically allowed region, 
which is delimited by the turning points, the classical probability density vanishes identically 
while the quantum-mechanical probability density is nonzero and decays exponentially. This 
discrepancy decreases as the energy increases because the quantum probability becomes more 
localized in the classically allowed region, but for any value of the energy there remains a 
sharp cutoff in the classical probability beyond the turning points. 

In the physical world the cutoff at the boundary between the classically allowed and the 
classically forbidden regions is not perfectly sharp. For example, in classical optics it is 
known that below the surface of an imperfect conductor, the electromagnetic fields do not 
vanish abruptly. Rather, they decay exponentially as functions of the penetration depth. 
This effect is known as skin depth [32] . The case of total internal reflection is similar: When 
the angle of incidence is less than a critical angle, there is a reflected wave and no transmitted 
wave. However, the electromagnetic field does cross the boundary; this field is attenuated 
exponentially in a few wavelengths beyond the interface. Although this field does not vanish 
in the classically forbidden region, there is no flux of energy; that is, the Poynting vector 
vanishes in the classically forbidden region beyond the interface. 

When classical mechanics is extended into the complex domain, classical particles are 
allowed to enter the classically forbidden region. However, in the forbidden region there is 
no particle flow parallel to the real axis. Rather, the flow of classical particles is orthogonal 
to the axis. This feature is analogous to the vanishing flux of energy in the case of total 
internal reflection as described above. 

We illustrate these properties of complex classical mechanics by using the classical har- 
monic oscillator, whose Hamiltonian is given in (jT|. The classical equations of motion are 

x = 2u, y = 2v, u = —2x, v = —2y, (2) 

where the complex coordinate is x + iy and the complex momentum is u + iv. For a particle 
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FIG. 3: Classical trajectories in the complex plane for the harmonic-oscillator Hamiltonian ([!]). 
These trajectories are nested ellipses. Observe that when the harmonic oscillator is extended into 
the complex domain, the classical particles may pass through the classically forbidden regions 
outside the turning points. When the trajectories cross the real axis, they are orthogonal to it. 



having real energy E and initial position x(0) = a> yE, y(0) = 0, the solution to (pq) is 



x(t) = acos(2t), y(t) = Va 2 -Esm(2t). (3) 

Thus, the possible classical trajectories are a family of ellipses parametrized by the initial 
position a: 

cr a 2 — E 

Four of these trajectories are shown in Fig. [3] Each trajectory has the same period T = ir. 
The degenerate ellipse, whose foci are the turning points at x = ±a/E, is the familiar real 
solution. Note that classical particles may visit the real axis in the classically forbidden 
regions \x\ > y/E, but that the elliptical trajectories are orthogonal rather than parallel to 
the real axis. 

We can use the solution in ^ to extend the plots of classical probabilities in Figs. [T] and 
[2] to the complex plane. The speed s of a classical particle at the complex coordinate point 
x + iy is given by 

s(x, y) = (E 2 + x 4 + y A + 2x 2 y 2 + 2y 2 E - 2x 2 E) 1/4 . (5) 

If we assume that all ellipses are equally likely trajectories, then the relative probability 
density of finding the classical particle at the point x+iy is l/s(x, y). This function is plotted 
in Fig. |4} Figure [4] shows that the classical-mechanical probability density extends beyond 
the classically allowed region and into the complex plane. In this figure the distribution of 
classical probability in the complex plane falls off as the reciprocal of the distance from the 
origin. Thus, beyond the turning points on the real line the classical probability density is 
no longer identically zero, and what is more, it begins to resemble the quantum-mechanical 
probability density in the classically forbidden regions on the real axis in Figs. [I] and |2j 

We emphasize that while the probability surface in Fig. [4] has the feature that it does 
not vanish on the real axis in the classically forbidden regions, it is entirely classical and 
thus there is none of the wavelike behavior that one would expect in quantum mechanics. 
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FIG. 4: Classical-mechanical relative probability density for finding a particle at a point x + iy in 
the complex plane for a particle subject to harmonic forces. The probability density resembles a 
pup tent with infinitely high tent poles located at the turning points. The complex classical ellipses 
in Fig. [3] are superposed on the tent canopy. 



This plot is two-dimensional and not one-dimensional, and this motivates us to define and 
calculate the complex analog of the quantum-mechanical probability density p(z) in the 
complex- ,2 plane. As in conventional quantum mechanics, this probability density p depends 
quadratically on the wave function ip, and by virtue of the time-dependent Schrodinger 
equation, p satisfies a local conservation law in the complex plane. 

Extending the quantum-mechanical probability density into the complex plane is non- 
trivial because i/j(z) is a complex-valued function of z, and as we will see, p(z) is not the 
absolute square of ip and thus p(z) is also a complex-valued function. In contrast with or- 
dinary quantum theory, it is not clear whether p(z) can be interpreted as a physical local 
probability density. To solve this problem we identify special curves in the complex- z plane 
on which p(z)dz is real and positive. We find these special curves by constructing a differ- 
ential equation that these curves obey. Discovering these curves then allows us to formulate 
the complex version of the correspondence principle. 

The complex correspondence principle is extremely delicate; understanding the prob- 
ability density in complex quantum mechanics requires sophisticated mathematical tools, 
including the use of asymptotic analysis beyond Poincare asymptotics; that is, asymptotic 
analysis of transcendentally small contributions. This paper is organized as follows: In Sec. [IT] 
we obtain the form of the complex probability density in complex VT quantum mechanics 
and discuss the specific case of the harmonic oscillator. Then, in Sec. Ill we give a simple 
example of a complex system, namely, a random walk with a complex bias, that has a real 
probability distribution in the complex plane. Next, in Sec. [TV] we study the special case of 
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the complex probability for the ground state of the quantum harmonic oscillator. In Sec. [V] 
we consider the complex probability for the excited states of the harmonic oscillator. The 



quasi-exactly-solvable "PT-symmetric anharmonic oscillator is discussed in Sec. VI Finally. 



in Sec. VII| we give some concluding remarks and discuss future directions for research. 



II. LOCAL CONSERVATION LAW AND PROBABILITY DENSITY IN THE 

COMPLEX DOMAIN 

In this section we derive a local conservation law associated with the time-dependent 
Schrodinger equation for a complex PT-symmetric Hamiltonian. Then we derive the formula 
for the local probability density and illustrate this derivation for the special case of the 
harmonic oscillator, which is the simplest Hamiltonian having VT symmetry. 

A. Local conservation law for VT quantum mechanics 

Consider the general quantum-mechanical Hamiltonian 

H = p 2 + V (x) . (6) 

The coordinate operator x and the momentum operator p satisfy the usual Heisenberg 
commutation relation [x,p] = i. The effect of the parity (space-reflection) operator V is 

VxV = —x and VpV = —p, (7) 

and the effect of the time-reversal operator T is 

TxT = x and TpT = —p. (8) 

Unlike V, which is linear, T is antilinear: 

TiT = -i. (9) 

Thus, requiring that H be VT symmetric gives the following condition on the potential: 

V*(-x) = V(x). (10) 

The time-dependent Schrodinger equation iijj t = Hip in complex coordinate space is 

ii> t {z, t) = t) + V(z)ij(z, t). (11) 



We treat the coordinate z as a complex variable, so the complex conjugate of (11) is 

-ir t {z%t) = -r*Az*>t)+v*(z*w(z*,t). (12) 

Substituting —z for z* in (fl2j), we obtain 



- irft(-z, t) = -i(>: s (-z, t) + V(z)r(-z, t), (13) 
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where we have used (10) to replace V*{— z) with V(z). To derive a local conservation law 



we first multiply (11) by z, t) and obtain 

iif>*(-z, t)il> t (z, t) = -if>*(-z, t)^ zz (z, t) + V(z)^*(-z, t)^(z, t). (14) 
Next, we multiply (FL3| by —ij)(z,t) and obtain 



ii>{z,t)r t {-z,t) = i>{z,t)r zz {-z,t) - V(zW(z,tW(-z,t). (15) 



We then add (14) to (15) and get 



- ty(z,t)i/>*{-z,t)) + — \ii/>(z,t)i/>l(-z,t) - ir(-z,t)M^t)} = 0. (16) 
This equation has the generic form of a local continuity equation 

p t (z,t) + j z (z,t) = (17) 

with local density 

p(z,t)=4>*(-z,t)if>(z,t) (18) 

and local current 

j( Z , t) EE ir z (~Z, t)i>{ Z , t) - iP(-Z, t)tl> g (z, t). (19) 



Note that the density p(z,t) in (18) is not the absolute square of ip(z,t). Rather, because 
VT symmetry is unbroken, if)*(—z,t) = if) z (z,t) and thus p(z,t) = [ip(z, t)] 2 . It is this fact 
that allows us to extend the density into the complex-z plane as an analytic function. 

B. Probability density for VT quantum mechanics 



While (17) has the form of a local conservation law, the local density p(z,t) in (18) is 
a complex-valued function. Thus, it is not clear whether p(z, t) can serve as a probability 
density because for a locally conserved quantity to be interpretable as a probability density, 
it must be real and positive and its spatial integral must be normalized to unity. With this 
in mind, we show how to identify a contour C in the complex- z plane on which p(z,t) can 
be interpreted as being a probability density. On such a contour p(z, t) must satisfy three 
conditions: 

Condition I : Im [p(z) dz] = 0, (20) 
Condition II : Re [p(z) dz] > 0, (21) 

Condition III : / p(z) dz = 1. (22) 
Jc 

A complex contour C that fulfills the above three requirements depends on the wave 
function ip(z,t) and thus it is time dependent. However, for simplicity, in this paper we 
restrict our attention to the wave functions ip(z,t) = e lEnt ip n (z) , where E n is an eigenvalue 
of the Hamiltonian and ip n {z) is the corresponding eigenfunction of the time-independent 
Schrodinger equation. For this choice of i/j(z,t) the local current j(z,t) vanishes, and p(z) 
and the contour C on which it is defined is time independent. (We postpone consideration 
of time-dependent contours to a future paper.) 
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This paper considers two PT-symmetric Hamiltonians, the quantum harmonic oscillator 

H = p 2 + x 2 (23) 

and the quasi-exactly-solvable quantum anharmonic oscillator [43] 

H = p 2 — x 4 + 2iax 3 + (a 2 - 2b)x 2 - 2i(ab - J)x. (24) 

For these Hamiltonians the eigenf unctions have a particularly simple form. The eigenfunc- 
tions of the harmonic oscillator are Gaussians multiplied by Hermite polynomials: 



- 2 / 2 He, 



z . 



(25) 



Analogously, the quasi-exactly-solvable eigenfunctions of the quartic oscillator in (24) are 



exponentials of a cubic multiplied by polynomials, as discussed in Sec. VI 



C. Quantum harmonic oscillator 

To find the probability density for the quantum harmonic oscillator in the complex-z 



plane, we use the eigenfunctions in (25) to construct p(z) according to (18) and then we 



impose the three conditions in (20) - (22). The function p(z) has the general form 



p(z,t) 



-x +y — 2ixy 



[S(x,y) +iT(x, 



where z = x + iy and S(x, y) and T(x, y) are polynomials in x and y: 

S(x,y) = Re ([He n (^)] 2 ) and T(x,y) = Im ([He^)] 2 ) . 
Thus, pdz has the form 

pdz = e~ x2+y2 [cos(2xy) — ism(2xy)][S(x,y) + iT(x, y)](dx + idy). 



(26) 



(27) 



(28) 



Imposing Condition I in (20), we get a nonlinear differential equation for the contour y(x) 
in the z = x + iy plane on which the imaginary part of p dz vanishes: 



dy S(x, y) sm{2xy) — T(x, y) cos(2x?/) 
dx S(x, y) cos(2xy) + T(x, y) sva(2xy) 



(29) 



On this contour we must then impose Conditions II and III in (21) and (22). To do so 
we calculate the real part of pdz: 



Re(pdz) = e x +v {[S(x, y) cos(2xy) + T(x, y) sm(2xy)\dx 
+ [S(x, y) sm.{2xy) — T(x, y) cos(2xy)]dy} . 



Thus, using (29), we get 



Re (p dz) = e 



[S(x,y)] 2 + [T(x,y)f 



S(x, y) cos(2xy) + T(x, y) s'm(2xy) 



dx 



(30) 



(31) 
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or, alternatively, 



Re (p dz) = e 



-x +y 



[S(x,y)} 2 + [T(x,y)f 



S(x, y) sm(2xy) — T(x, y) cos(2xy) 



dy. 



(32) 



Equations (31) and (32) suggest that there may be some potentially serious problems 



with establishing the existence of a contour y(x) in the complex plane on which p could be 
interpreted as a probability density. First, it appears that if the contour y(x) should pass 
through a zero of either the numerator or the denominator of the right side of ( 29 ) , then 



the sign of Re (pdz) will change, which would violate the requirement of positivity in (21) 



However, in our detailed mathematical analysis of these equations in Sec. [TV] we establish the 
surprising result that the sign of Re (p dz) actually does not change. Second, if the contour 
y(x) passes through a zero of the numerator or the denominator of the right side of (29), one 
might expect p to be singular at this point. Indeed, it is singular, and one might therefore 



worry that the integral in (22) would diverge and thus violate condition III that the total 



probability be normalized to unity. However, we show in Sec. IV that the singularity in p is 
an integrable singularity, and thus the probability is normalizable. 



III. COMPLEX RANDOM WALKS 

To illustrate at an elementary level the concept of a local positive probability density along 
a contour in the complex plane, we consider in this section the case of a one-dimensional 
classical random walk determined by tossing a coin having an imaginary bias. For a conven- 
tional one-dimensional random walk, the walker visits sites on the real-x axis, and there is a 
probability density of finding the random walker at any given site. However, for the present 
problem the random walker can be thought of as moving along a contour in the complex-z 
plane. For this introductory problem the contour is merely a straight line parallel to the 
real-x axis. 

In the conventional formulation of a one-dimensional random-walk problem the random 
walker may visit the sites x n = nS, where n = 0, ±1, ±2, ±3, ... on the x axis. These sites 
are separated by the characteristic length 5. At each time step the walker tosses a fair coin 
(a coin with an equal probability 1/2 of getting heads or tails). If the result is heads, the 
walker then moves one step to the right, and if the result is tails, the walker then moves one 
step to the left. To describe such a random walk we introduce the probability distribution 
P(n, k), which represents the probability of finding the random walker at position x n at time 
tfc = kr, where k = 0, 1, 2, 3, . . .. The interval between these time steps is the characteristic 
time t. The probability distribution P(n, k) satisfies the partial difference equation 

Pn,k — \Pn+l,k-l + \Pn-l,k-l- (33) 

Let us generalize this random-walk problem by considering the possibility that the coin 
is not fair; we will suppose that the coin has an imaginary bias. Let us assume that at each 
step the "probability" of the random walker moving to the right is 

a = \ + i@8 (34) 

and that the "probability" of moving to the left is 



1 - a = | — i/35 = a*. 



(35) 
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The parameter f3 is a measure of the imaginary unfairness of the coin. Note that a + a* = 1, 
so the total "probability" of a move is still unity. Now, (J33|) generalizes to 



P(n, k) = aP(n - 1, k - 1) + a*P{n + 1, k - 1). (36) 

We wish to interpret P(n, k) as a probability density, but such an interpretation is un- 
usual because P{n, k) obeys the complex equation (36) and is consequently complex-valued. 
Although this equation is complex, it is VT symmetric because it is invariant under the 
combined operations of parity and time reversal. Parity interchanges the probabilities (34) 
and (35) for rightward and leftward steps and thus has the effect of interchanging the indices 
n+ 1 and n — 1 in (36). Time reversal is realized by complex conjugation. 

We choose the simple initial condition that the random walker stands at the origin n = 
at time k = 0: 

P(0,0) = 1. (37) 



For this initial condition the exact solution to <\36\j is 

P(n, k) 



k\oT^ (a*)*" 2 " 



(*±») !(*=*)! 

(38) 



, . k+n , . k — n 

kla 2 l — a) 2 



This solution is manifestly VT symmetric because it is invariant under combined space 
reflection n — > —n and time reversal (complex conjugation). 

Now let us find the continuum limit of this complex random walk. To do so we introduce 
the continuous variables t and x according to 

t = kr and x = n5 (39) 

subject to the requirement that the diffusion constant a given by the ratio 

a = 5 2 /t (40) 

is held fixed. We then let S — > and define p(x, t) to be the probability density: 

P(n,k) 1 _ (*-2»<r/3t) 2 

p( ,, () , te _ = _ e « . (41) 

The function p(x, t) solves the complex diffusion equation 

p t (x, t) = o-p xx (x, t) - 2io-(3p x (x, t) (42) 

and satisfies the delta-function initial condition lim^o p{x,t) = S(x), where x is real. 

Observe that there is a contour in the complex-x plane on which the probability density 
p(x, t) is real. For this simple problem the contour is merely the straight horizontal line 

lmx = 2o-/3t. (43) 



At t = this line lies on the real axis, but as time increases, this line moves vertically at the 
constant velocity 2o~(3. Thus, even though the probability density is a complex function of 
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x and t, we can identify a contour in the complex- x plane on which the probability is real 
and positive and hence may be interpreted as a conventional probability density. 

The situation for complex quantum mechanics is not so simple. For the quantum problem 
we seek a complex contour that satisfies the three conditions (20) - (22). On such a contour 
the probability density p(z, t) for finding a particle in the complex-z plane at time t is in 
general not real. Rather, it is the product p dz (which represents the local contribution to the 
total probability) that is real and positive. In the case of the complex random-walk problem 
in this section the contour happens to be a horizontal line, and thus both the infinitesimal 
line segment dx and the probability density p are individually real and positive. 



IV. PROBABILITY DENSITY ASSOCIATED WITH THE GROUND STATE OF 

THE HARMONIC OSCILLATOR 

In Sec. [IT] we showed that for the quantum harmonic oscillator, the differential equation 



(29) defines a complex contour C on which pdz, the local contribution to the total proba- 



bility, is real and positive. In this section we discuss the mathematical techniques needed to 



analyze differential equations of this form. Surprisingly, even the simplest version of (29), 

dy 

— = sin(2xy)/ cos(2xy), (44) 

(JjJb 

which is associated with the ground state of the harmonic oscillator, does not have a closed- 



form solution. Note that (44) is the special case of (29) for which S = 1 and T = 0. 

First-order differential equations of the general form y'(x) = f(ax + /3y) (where a and f3 
are constants) and the general form y'(x) = f(y/x) have simple quadrature solutions, but 
finding exact closed-form solutions to differential equations of the general form y'(x) = f(xy) 



is hopeless. Analyzing the behavior of solutions to equations such as (44) requires the use 
of powerful asymptotic techniques explained in this section. 



A. Toy model 

We begin by examining a simplified toy- model version of (44), namely, 

^ = cos(xy), (45) 

which is given as a practice problem in Ref. jS]. Numerical solutions to this equation 
for various initial conditions y(0) = 1.0, 2.0, 3.0, 4.0, 4.5, 5.0, 5.5, 6.0 are shown in Fig. [5j 
Observe that the solutions become quantized; that is, after leaving the y axis, they bunch 
together into isolated discrete strands, which then decay smoothly to 0. The first bunch of 
solutions all have one maximum, the second bunch all have two maxima, and so on. 



We cannot solve (45) exactly, but we can use asymptotics to determine the features of the 
solutions in Fig. [5] To do so, we must first examine the behavior of the quantized strands 
for large x. Note that each strand becomes horizontal as x — > oo and that the solutions to 
the differential equation have zero slope on the hyperbolas 

xy = (n + |) 7T, (46) 
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yix) 




FIG. 5: Solutions to the differential equation y'(x) = cos(xy) for various initial conditions. Note 
that the solutions bunch together into discrete strands as x increases. 



where n is an integer (see Fig. [6]). This suggests that the possible leading asymptotic behavior 
of y(x) for large x is given by 



y(x) 



(» + §) 



7T 



[X — > OO 



X 



(47) 



To verify (47) one must determine the higher-order corrections to this asymptotic behavior. 



These corrections take the form of a series in inverse powers of x: 



y{x) 



(n + |) 7T 



X 



+( 



1 3 , N 90+ (n + |V 

n _|_ _|_ I j\n \ 2/ 



1 + (-1)" 1 h (-1)" — V 



7T 2 315 + 8 (n + |) 



2 2 
7T 



X 2 X 4 



37800 + 1440 (n + l) 2 TT 2 + 3(n+l)' 
' ' 40x 10 



7T 



+ 



3x 8 



(48) 



A remarkable feature of the asymptotic behavior in (48) is that there is no arbitrary 



constant. In general, the exact solution to an ordinary differential equation of order iV must 
contain iV arbitrary constants of integration and the asymptotic behavior of the solution 
must also have iV arbitrary constants. (The constant n is not arbitrary because it is required 
to take on integer values.) Thus, the paradox is that the complete asymptotic description of 



the solution to (45) must contain an arbitrary constant, and yet no such constant appears 



to any order in the asymptotic series in (48) 



The resolution of this puzzle is that there is an arbitrary constant of integration in the 
asymptotic behavior (48), but it cannot be seen at the level of Poincare asymptotics. In 
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yixi 




2 4 6 8 10 



FIG. 6: Solutions (solid lines) to y'(x) = cos(xy) for various initial conditions. The slope of y(x) 
vanishes on the hyperbolas (dotted lines) xy = (n + ^) ir, where n = 0, 1, 2, .... Asymptotic 
analysis beyond all orders explains why solutions bunch together and approach the even-numbered 
hyperbolas. When n is odd, the curves near the hyperbolas are unstable and veer away from the 
hyperbolas. To illustrate this instability, three initial conditions on either side of the separatrix 
initial condition y(0) = 5.276 032 283 736 901518 295442... are shown. The solutions beginning 
at y(0) = 5.287 and at y(0) = 5.266 veer off between x = 4 and x = 5, those beginning at 
y(0) = 5.276 032 296 and at y(0) = 5.276 032 267 veer off near x = 7, and those beginning at 
y(0) = 5.276 032 283 736 901 784 and at y(0) = 5.276 032 283 736 901 484 veer off near x = 9. 



Poincare asymptotics one ignores transcendentally small (exponentially small) contributions 
to the asymptotic behavior because such contributions are subdominant (negligible compared 



with every term in the asymptotic series). The missing constant of integration in (48) can 
only be found at the level of hyperasymptotics (asymptotics beyond all orders) [4"5]44"T] . 
To find the missing constant of integration, we first observe that the difference of two 



solutions to (45) that belong to different bunches (that is, corresponding to different values 
of n) approaches zero like 1/x for large x: 



ym{x) - y n (x) 



(m — n)7t 



[X — > oo 



X 



(49) 



However, the difference between two solutions belonging to the same bunch is exponentially 
small. To see this, let y(x) and u(x) be two solutions in the nth bunch and let D(x) = 
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y(x) — u(x) be their difference. Then D(x) obeys the differential equation 



dD 

dx 



= cos(xy) — cos(xu) 

= — 2 sin [|:r.D(a;)] sin [^x(y + u) 

xD(x) sin [(n + I) vr + (-1) 71 (n + |) tt/x 2 + . . .] (x -»> oo) 

~ -x£>(x)(-l) n cos [(n + |) 7r/x 2 + . . .] (x -> oo) 

arD(a;)(-l) n [l - § (n + ±) 2 7r 2 /a; 4 + ■ ■ -1 -> oo), 



(50) 



where we have assumed that D(x) is exponentially small and hence that the average of y(x) 
and u{x) is given by (48). Thus, the solution for D{x) is 



D(x) ~ Ce-^ 1 )"* 3 / 2 



1 - o 

v ; Ax 2 



+ 



(x oo), 



(51) 



where C is an arbitrary multiplicative constant of integration because the differential equa- 



tion (50) is linear. 



The derivation of (51) is valid only when n is an even integer because the function D(x) 



is exponentially small only when n is even. When n is odd, the result in (51) is of course 



not valid; however, we learn from this asymptotic analysis that the asymptotic solutions in 



(48) are unstable for odd n. That is, nearby solutions veer off as x increases and become 
part of the adjacent bunches of solutions (see Fig. [6]). One solution, called a separatrix, lies 
at the boundary between solutions that veer upward and solutions that veer downward; the 



separatrices are the solutions whose asymptotic behaviors are given by (47) and (48) for odd 
n. On Fig. [6] the separatrix corresponding to n = 5 is shown; this separatrix evolves from 
the initial condition y(0) = 5.276 032 283 736 901 518 295 442 . . .. 



B. Stokes' wedges and complex probability contours 

The complex probability density p(z) for the quantum harmonic oscillator is the square 
of the eigenfunctions of the time-independent Schrodinger equation. These eigenfunctions 
vanish like e~ z I 2 as z — > ±oo on the real axis. More generally, these eigenfunctions vanish 
exponentially in two Stokes' wedges of angular opening n/2 centered on the positive-real and 
the negative-real axes. We refer to these Stokes' wedges as the good Stokes' wedges because 



this is where the probability integral (22) converges. Conversely, the eigenfunctions grow 
exponentially as \z\ — > oo in two bad Stokes' wedges of angular opening 7r/2 centered on the 
positive-imaginary and negative-imaginary axes. The probability integral clearly diverges if 
the integration contour terminates in a bad wedge. 

Our objective is to find a path in the complex- ,2 plane on which the probability is real 



and positive. This path then serves as the integration contour for the integral in (22). The 
total probability can only be normalized to unity if this integral converges. Therefore, this 
contour of integration must originate inside one good Stokes' wedge and terminate inside 
the other good Stokes' wedge. 

When we began working on this problem, we were surprised and dismayed to find that, 
apart from the trivial contour lying along the real axis, a continuous path originating and 
terminating in the good Stokes' wedges simply does not exist! The (wrong) conclusion that 
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one might draw from the nonexistence of a continuous integration path is that it is impossible 
to extend the quantum-mechanical harmonic oscillator (and quantum mechanics in general) 
into the complex domain without losing the possibility of having a local probability density. 
To investigate the provenance of these difficulties it is necessary to perform an asymptotic 



analysis of the differential equation (44) that defines the integration contour. Using the 
techniques explained in the previous subsection, we find that deep inside the Stokes' wedges, 
the solutions to this differential equation approach the centers of the wedges. Specifically, 
in the good Stokes' wedge 

T17T 

V(x) ~ — (x ->■ +00), (52) 
where n is an integer, while inside the bad Stokes' wedge 

x{y)~ { ^^- (y^+00), (53) 

where m is an integer. Thus, as \z\ — > 00, the integration paths become quantized in the 
same manner as solutions to the toy-model discussed above. However, solutions deep in the 
good wedges are all unstable (like the odd-n solutions in the toy model), while solutions 
deep in the bad wedges are all stable (like the even-n solutions in the toy model). Thus, as 
\z\ — > 00 in a bad wedge, all solutions combine into quantized bunches and then approach 



the curves in (53). On the other hand, as \z\ — > 00 in a good wedge, all solutions except for 
the isolated separatrices in (52) peel off, turn around, and leave the Stokes' wedge. These 
solutions then approach infinity in one of the bad Stokes' wedges. This behavior is illustrated 
in Fig. [7j 



Because of the inherent instability in the differential equation (44), only a discrete set of 
separatrix solutions actually enter and remain inside the good Stokes' wedges. Finding these 
special curves is the analog of calculating the eigenvalues of a time-independent Schrodinger 
equation. It is only possible to satisfy the boundary conditions on the eigenfunctions of 
a Schrodinger equation for a discrete set of eigenvalues. Thus, by analogy we could refer 
to these special separatrix paths in the complex plane as eigenpaths jlH]. With the single 
exception of the path that runs along the real axis, all eigenpaths terminating in a good 
Stokes' wedge must originate in a bad Stokes' wedge, as shown in Fig. [7j 

Thus, it appears that apart from the real axis there is no path along which the probability 



integral (22) converges. Despite this apparently insurmountable problem, it is actually 



possible to find contours in the complex- ,2 plane on which the probability density satisfies 



the three requirements (20) - (22), as we explain in the next subsection. 



C. Resolution of the problem posed in Subsection IV B 



Amazingly, it is possible to overcome the difficulties described in Subsec. |IVB In this 
subsection we use asymptotic analysis beyond all orders to show that if the probability 



integral (22) is taken on a contour C\ that enters a bad Stokes' wedge and continues off to 



00 in the wedge and then the integration path emerges from the bad wedge along a different 
contour C2, the combined integral along C\ + C2 exists if C\ and C2 are solutions to (44). 



This conclusion is valid even though the integrals along C\ and along C2 are separately very 
strongly (exponentially^ divergent. This is the principal result of this paper. 



Recall that on every contour that solves the differential equation (44), the infinitesimal 



probability p(z) dz is real. Now the objective is to select from all these contours those special 
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FIG. 7: Numerical solutions (solid curves) to the differential equation (44) in the complex- z = x+iy 
plane. These solutions have vanishing slope on the lightly dotted hyperbolas and infinite slope on 
the heavily dotted hyperbolas. There are two "good" Stokes' wedges (dark shading) inside of 
which p(z) decays to exponentially as \z\ — > oo along the solution curves and there are two "bad" 
Stokes' wedges (unshaded) inside of which p(z) grows exponentially as \z\ — > oo along the solution 
curves. Note that as a typical solution curve goes more deeply into a good Stokes' wedge, the 
curve becomes unstable and eventually turns around. The solution curve then leaves the good 
Stokes' wedge and is pulled around into a bad Stokes' wedge. In the bad Stokes' wedge curves are 
stable and continue on to ±ioo. The only curves that ever reach ±oo in the good Stokes' wedges 
are separatrices. Five such separatrix paths labeled n = 0, n = 1, ...,n = 4 are shown as heavy 
dashed curves. One special separatrix path runs along the real axis and connects the good Stokes' 
wedge at — oo to the good Stokes' wedge at +oo. This is the only continuous unbroken curve that 
connects the two good Stokes' wedges. Eight bunches of curves (labeled by m) that enter the bad 
Stokes' wedges are shown. 



contours on which the integral (22) exists. For the ground state of the harmonic oscillator 
this integral becomes 

dye v\-[<y)? 1 (54) 
sm[2yx[y)\ 

To analyze the integral (54) we must use the Poincare asymptotic series expansion of x(y) 
for large y, 



(™±ik d + J_ + A + 45 ~( m + 1) 

2y I 2y 2 4y 4 24y 



o. l 1 + ^ + T^ + -^T 1 + •••) (v^+oo), (55) 
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whose leading asymptotic behavior is given in (53). The coefficients in this asymptotic 
series are easily derived from the differential equation (44). The asymptotic behavior of 
x(y) beyond all orders is also required. To find this behavior, we let u(y) and v(y) be two 
solutions to ^ = cos(2xy)/ sm(2xy) having the same value of n (that is, belonging to the 
same bunch). The Poincare asymptotic behavior of these solutions is given in (55), but if 
we define D(y) = u(y) — v(y), we find that 



D(y) ~ Ce- y 



(m+\)' 



7T 2 12(m+i)V + (m + f)V 



Ay 2 



+ 



32y* 



+ 



(y +oo) 



(56) 

where C is an arbitrary constant. 

Consider an integration contour that enters the bad Stokes' wedge centered on the positive 
imaginary axis. This contour starts at y = Y, follows the path x = u(y), and runs up to 
y = +oo. The contour then leaves the wedge along the path v(y), and it returns to y = Y. 
Such a contour is shown in Fig. [8j The total contribution along both paths u(y) and v(y) 
to this integral is 



Y 



Y 



dy e yZ 
dy e y * 



Hy)] 2 



My)? \ 



sva[2yu{y)\ sin[2yv(y)] J 

sm[2yv(y)}e- [u(y)]2 - sm[2yu(y)]e-M y)]2 
sm[2yu(y)] sm[2yv(y)) 



(57) 



Let us now investigate the convergence of this integral. We approximate the denominator 
of the integrand by using ordinary Poincare asymptotics in which we ignore the transcen- 
dentally small difference between u(y) and v(y). We first establish that 



2yu{y) ~ 2yv{y) 



sm 



. (m+ |) 7r 3 (m + I) tt 

(m + \) tx + v - r + 



2y 2 



Ay* 



cos 



(m + |) 7T 3 (m + |) 7T 



2y* 

i m + \)'' 
8y 4 



Thus, the denominator is given approximately by 



+ 



7T 



Ay* 



+ ... 



+ 



(y -»■ +oo) 



(58) 



sm[2yu(y)]sm[2yv(y)} 



Ay* 



7l~ 



o 



(y -»■ +oo) 



(59) 



Next, we approximate the numerator of (57). Since u(y) and v(y) are small as y — > +oo, 



we expand the exponentials. Keeping two terms in the expansion, we find that 

sm[2yv{y)]e- Hy)]2 - sm[2yu{y)]e- [v(y)]2 

~ sm[2yv(y)} (l - [u{y)] 2 + ...)- sm[2yu(y)} (l - [v(y)} 2 + . . .) 



(60) 
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FIG. 8: A contour in the complex- z plane comprised of two numerical solutions to (44). The first 
path starts on the positive-imaginary axis at y = 0.003 324973 872 707 912 (the starting point is 
indicated by a dot), and after it becomes vertical, it veers off to y = oo in the bad Stokes' wedge. 
The second path is the separatrix belonging to the same bunch as the first path. This separatrix 
leaves the bad Stokes' wedge and runs to x = oo in the good Stokes' wedge. Even though the 



probability density blows up exponentially in the bad Stokes' wedge, it is shown in (65) that the 
integral of the probability density from the value y = Y up to y = oo along the first path and then 
from y = oo down to y = Y along the second (separatrix) path is convergent. 



as y — > +oo. Let us see what happens if we neglect the terms of order u 2 and v 2 and keep 



only the first pair of terms in (60). Using the trigonometric identity 



sin a — sin (3 = 2 cos 



a + (3 



sm 



a — (3 



(61) 
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we find that the numerator has the following hyperasymptotic form: 
sm[2yv(y)] — sm[2yu(y)] ~ 2 cos 



, , , (m + I) 7r 3 (m + |) 
(m+ §)tt + v - 27 + 



7T 



2y 2 



V 



sin(— yD) 



2yD(y)(-l) m sm 



(m + \) n 3 (m + §) vr 



2^ 



(-1) 7 



y 



, 3 (m + |) 7r / 1 

(m + I) 7r + v _ „ 2/ + O 



(62) 



Hence, the leading-order contribution to the integrand of the integral in (57) gives 

Mv) 



d ye y (-iy 
y y 



(m + |)7r. 



(63) 



Thus, apart from an overall multiplicative constant, this integral has the form J^dy/y. 
This integral is not exponentially divergent, but it is still logarithmically divergent and thus 
it is not acceptable. 



Fortunately, this logarithmic divergence cancels if we include higher-order terms in (60) 



in the calculation; including the quadratic terms in the expansion of the exponentials is 



sufficient to make the integral (57) converge. In addition to the terms that we examined in 



(62), we also consider 



[v(y)] 2 sm[2yu(y)} - [u{y)f sva[2yv{y)\ 
= ({v(y)] 2 - [u(y)] 2 ) sm[2yu(y)] + [u 

'1 1 



-n(-l) m (m+±)D(y) 
+ (-l) m (m+!)Vl>(y) 



y 2 y 3 
l 

% 3 



o 



2 (sxn[2yu{y)) - em[2yv(y)]) 
1 



+ 



(y +oo). 



(64) 



We combine this asymptotic contribution to the numerator of the integrand of (57) with the 



leading-order contribution in (62) and obtain a further cancellation. The resulting integral 



dy (-l) m e^ 



D(y) 



y 3 



71 



C(-l) m (m+i) 



y y 



dy 

3 



(65) 



converges! Because this integral is convergent, a contour can begin at a point on the positive 
imaginary axis, run into and back out of the bad Stokes' wedge at y = +oo, and then 
terminate in the good Stokes' wedge at x — +oo. Such a contour, which begins at y = 
0.003 324973 872 707912 on the imaginary axis, is shown in Fig. [9} 

Finally, if we reflect this contour about the y axis, we get a complete path originating 
in the left good Stokes' wedge at x — — oo, crossing the y axis horizontally, and eventually 
terminating in the right good Stokes' wedge at x — +oo. Such a contour is shown in the 
upper plot in Fig. [9j Note that a contour that leaves the bad Stokes' wedge need not leave 
on a separatrix path that terminates in the good Stokes' wedge; instead, it can leave the bad 
Stokes' wedge and follow a path that turns around and returns to the bad Stokes' wedge. 
Such a path, which is shown in the lower plot in Fig. [9j finally leaves the bad Stokes' wedge 
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FIG. 9: Complex contours that run from one good Stokes' wedge to the other. In the upper plot a 
contour begins at x = — oo in the left good Stokes' wedge (shaded region), leaves the good Stokes' 
wedge along a separatrix path, and runs up to zoo in the bad Stokes' wedge (unshaded region). 
The contour then continues downward along a path in the same bunch, crosses the imaginary axis 
at y = 0.003 324 973 872 707912, and heads upwards into the same bad Stokes' wedge. Finally, the 
contour re-emerges from the bad Stokes' wedge and continues towards x = +oo along a separatrix 
path in the right good Stokes' wedge. Contours have zero slope on the lightly dotted lines and 
infinite slope on the heavily dotted lines. In the lower plot the contour connecting the two good 
Stokes' wedges visits the bad Stokes' wedge four times instead of twice. 

on a separatrix path and terminates in the good Stokes' wedge. Thus, the contour can be 
multiply thatched before it eventually terminates in the good Stokes' wedge. 

We remark that the device described here, in which the contour enters and then leaves a 
bad Stokes' wedge, does not work with the method of steepest descents, which is a standard 
technique used to find the asymptotic behavior of integrals. When a steepest path runs off 
to oo, it can only do so in a good Stokes' wedge and not in a bad Stokes' wedge; otherwise, 
the integral would diverge. In contrast, for the problem discussed in this paper the path 
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of real probability has no choice; it necessarily enters the bad Stokes' wedge, and it is the 
stable bunching of contours into quantized strands that saves the day. 



D. Safely crossing lightly dotted lines and heavily dotted lines 



There is one more potential problem that must be considered before we can claim to 
have a complete contour on which the probability density is positive and normalizable. It 
is necessary to show that the sign of the probability density does not change if the contour 
crosses one of the hyperbolas 2xy = nn on which y'(x) = 0, or 2xy = (m + |) n, on which 
y'(x) = oo. These hyperbolas are shown as lightly dotted and heavily dotted lines in Fig. [9j 

We can see that the contours in Fig. [9] cross both lightly and heavily dotted lines. When 
this happens, the sign of the denominators in (31) or (32) change, so at first one might 



think that the probability density along the contour would not remain positive. To address 



this concern, we rewrite each of the differential probabilities (31) and (32) as a differential 



probability proportional to the infinitesimal path length element ds by using the standard 
formula ds = a/ dx 2 + dy 2 . Thus, in (31), for example, we eliminate dx in favor of ds and 
obtain 

(66) 



dx 



We then substitute the differential equation (29) into (66) and get 



Re(pdz) = e~ x2+y2 y/\S{x, y)} 2 + [T(x, y)} 2 ds. 



(67) 



Clearly, this form of the infinitesimal probability contribution along the contour is explicitly 
positive and cannot change sign. 

Where is the error in the reasoning that led us to worry that the sign of Tle(pdz) might 
change? We note that the contour in Fig. [9] crosses a heavily dotted line near x = 2 and 
y — 1/2. Thus, the denominator of (31 ) does indeed change sign. However, as the probability 



curve becomes vertical, it simultaneously changes direction (that is, it runs backward). This 
is equivalent to dx changing sign, and this change in sign compensates for the change in sign 
of the denominator. Thus, all along the contour the infinitesimal contributions to the total 
probability remain positive. 



V. COMPLEX PROBABILITY DENSITY FOR HIGHER-ENERGY STATES OF 

THE HARMONIC OSCILLATOR 



In this section we generalize the analysis of Sec. IV to the excited states of the quantum 
harmonic oscillator. The eigenfunctions of these states differ from the ground state in that 



they have nodes on the real axis. As a consequence, the functions S(x, y) and T(x, y) in (26 ) 
become more complicated and as a result the differential equation (29), which determines 



the probability eigenpaths, is correspondingly more challenging to analyze. 



The first excited state ip\( 



ze 



2 12 



whose energy is 3, has a single node, which is 



located at the origin. For this eigenfunction S(x, y) — x — y and T(x, y) = 2xy. Thus, the 



differential equation (29) becomes 



dy 

dx 



X' 



— y 2 ) sin(2x?/) — 2xy cos(2xy) 
[x 2 — y 2 ) cos(2xy) + 2xysm{2xy) 



(68) 
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It is necessary to determine the asymptotic behavior of solutions to this equation near the 
node at the origin in the (x,y) plane. To do so we seek a leading asymptotic behavior of 
the form y(x) ~ ax + . . . as x — > 0. Substituting this behavior into (68) gives an algebraic 
equation for a: 

2a , , 

a = (69) 

The solutions to this equation are a = and a = ±a/3, which indicates that eigenpaths 
can enter or leave the node in one of six possible directions separated by 60°. (These paths 
are indicated on Fig. 10, except that the trivial path along the real axis has been omitted.) 
Apart from the behavior in the vicinity of the node, the eigenpaths shown in this figure are 
qualitatively similar to those shown in Figs. [7] and [9j 

shows that there are many compound eigenpaths connecting the two good 



Figure 10 



Stokes' wedges in addition to the conventional path along the real axis. A typical eigenpath 
begins in the left good Stokes' wedge and runs along a separatrix into one of the bad Stokes' 
wedges centered on the positive- or negative- imaginary axes. The path then emerges from 
and returns to the bad Stokes' wedge several times before crossing the imaginary axis. This 
path crosses the imaginary axis in two possible ways: (i) The path may cross the imaginary 
axis at the node, and if it does, the path may form a cusp, (ii) The path may cross the 
imaginary axis at a point other that is not a node, in which case it must be horizontal at the 
crossing point. The path then enters and emerges from a bad Stokes' wedge several more 
times before running off along a separatrix to infinity in the right good Stokes' wedge. 

The eigenfunction representing the second excited state has the form ifriiz) = {2z 2 — 
l)e~ z I 2 . The energy is 5. There are now two nodes, which are located on the real axis 
at ±l/\/z. These nodes are shown on Fig. 11 The eigenpaths in Fig. |ll| are qualitatively 



similar to those shown in Fig. 10 Like the eigenpaths shown in Fig. [LOl the eigenpaths in 



Fig. 11 enter the nodes horizontally or at ±60° angles. 



VI. COMPLEX PROBABILITY DENSITY FOR THE QUASI-EXACTLY 
SOLVABLE ANHARMONIC OSCILLATOR 

In this section we generalize the results of the previous two sections for the quantum 
harmonic oscillator to the more elaborate case of the quasi-exactly-solvable PT-symmetric 



anharmonic oscillator, whose Hamiltonian is given in (24). We consider here the special 
class of these Hamiltonians for which 6 = because the case b ^ presents no additional 
features of interest. The time- independent Schrodinger equation for this Hamiltonian is 

d 2 \ 

x 4 + 2iax 3 + a 2 x 2 — 2iJx \ if) n (x) = E n if) n (x). (70) 



dx 2 

When J is a positive integer, the first J eigenfunctions have the form 

^(x) = e- ix3 / 3 - ax2 / 2 P J - 1 (x), (71) 

where 

J-2 

P J _ 1 (x)=x J - 1 + ^2c k x k (72) 



fc=0 

is a polynomial of degree J — 1. 
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FIG. 10: Complex contours that run from the left good Stokes' wedge to the right good Stokes' 
wedge for the first excited state of the quantum harmonic oscillator. Four contours (solid lines) 
that begin at x = — oo in the left good Stokes' wedge (shaded region) are shown. These contours 
leave this Stokes' wedge along separatrix paths and run off to ±ioo in the upper and lower bad 
Stokes' wedges (unshaded regions). After visiting a bad Stokes' wedge one or more times, the 
contours pass through the node at the origin at 60° angles to the horizontal. At this node the 
probability density vanishes. Then the contours repeat the process in the right-half plane and 
eventually enter the right good Stokes' wedge along separatrix paths. Note that it is also possible 
to have a complex contour that does not pass through the node at the origin and still connects the 
good Stokes' wedges. The solution curves are horizontal on the lightly dotted lines and vertical on 
the heavily dotted lines. 



A. Probability density in the complex plane for the ground state 

We begin the analysis by considering the case J = 1 for which the ground-state eigen- 
function can be found exactly and in closed form: 

a -ix 3 /3-ax 2 /2 



(73) 



The associated ground-state energy is E = a. According to (|18|), the time-independent 

p( z ) = e -a*7»-«*-. (74) 



local probability density p in the complex-;? plane for this eigenfunction is 

-2iz 3 /3-az 2 
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FIG. 11: Complex contours that run from one good Stokes' wedge to the other for the case of the 
second excited state of the quantum harmonic oscillator. Two contours (solid lines) that begin 
at x = — oo in the left good Stokes' wedge (shaded region) are shown. These contours leave the 
left good Stokes' wedge along separatrix paths and run off to ±ioo in the upper and lower bad 
Stokes' wedges (unshaded regions). After visiting the bad Stokes' wedge one or more times, the 
contours pass through the left node on the real axis at 60° angles to the horizontal. At this node 
the probability density vanishes. The contours then reenter the upper or lower bad Stokes' wedges 
and pass through the right node on the real axis. After returning to the bad Stokes' wedges yet 
again, the solution curves eventually enter the right good Stokes' wedge along separatrix paths. 
The solution curves are horizontal on the lightly dotted lines and vertical on the heavily dotted 
lines. 



As in our study of the harmonic oscillator in Sees. IV and |Vj we take z = x + iy and 
dz = dx + idy. We then obtain 

».2„, o„.3 /r>_ r ,~.2 i „„. 2 



p{ Z )dz = e 2x W-V/3-«* +ay _ isin Q^ dx + ^/), 



(75) 



where 



9 = -2xy 2 + 2x 3 /3 + 2axy. 



(76) 

Consequently, Condition I in (20), which requires that Im [p(z)dz] = 0, translates into the 
differential equation 

dy = sing 
dx cos 9 
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FIG. 12: Good Stokes' wedges (shaded regions) and bad Stokes' wedges (unshaded regions) for 
the quasi-exactly-solvable VT anharmonic oscillator in (70). In the three good Stokes' wedges the 
eigenfunctions decay to zero like the exponential of a cubic [see (73)]. The eigenfunctions grow like 
the exponential of a cubic in the bad Stokes' wedges. 



This differential equation is the analog of (44) for the harmonic oscillator. Also, Condition 



III in (22), which requires that the integral of Re [p(z)dz] exist, is the same as demanding 
that the following (equivalent) integrals exist: 



2x 2 y— 2j/ 3 /3— ax 2 +ay 2 



2x 2 y-2y 3 /3-ax 2 +ay 2 



dy_ 
sin 9 



d X 2x 2 y-2y 3 /3-ax 2 +ay 2 



COS 9 



(78) 
(79) 



The integration contours for the integrals above must terminate in good Stokes' wedges 
in order that the integrals converge. The locations and opening angles of the Stokes' wedges 
are identified by determining where the probability density p(z) is exponentially growing 
or decaying. There are six Stokes' wedges, each having angular opening n/3. Of the three 
good Stokes' wedges [where p(z) decays exponentially] one is centered about the positive- 
imaginary axis and the other two lie adjacent to and below the positive-real and negative- real 
axes. One of the three bad Stokes' wedges is centered about the negative- imaginary axis 
and the other two lie adjacent to and above the positive-real and negative-real axes. The 



good and bad Stokes' wedges are shown in Fig. 12 as shaded and unshaded regions. 



1. Asymptotic analysis of solutions in the good Stokes' wedge below the positive-real axis. 

We begin by finding the asymptotic behavior of a function y(x) that solves the differential 
equation (77) and which approaches the center of the good wedge that lies adjacent to 
and below the positive-real z axis. (Because of VT symmetry, which is simply left-right 
symmetry in the complex-z plane, the left good wedge is treated in a similar fashion.) The 
center of this wedge lies at an angle of —tt/6 in the complex-;? plane, so we seek a solution 



26 



y(x) that satisfies the asymptotic condition 

y ~ -x/V3 (x -> +00). (80) 
The full asymptotic series approximation for such a solution has the form 

x A 2 A 3 Ai A 5 A 6 

y — ^ + A i + — + ^ + ^3 + ^ + ^ + ---- ( 81 ) 

To determine the coefficients in this series we observe that if we choose 

A 1 = \a and A 2 = -\a 2 V?>, (82) 



then 9 -»■ 4A 3 /\/3 as a; —7- 00. Then, using the differential equation rt77J ), we get —l/y/l = 
tan(4A3/v3), which gives 

A 3 = 4(n-l)n. (83) 



4 V 6 

With this choice of A 3 , we can use the differential equation to determine all of the higher- 
order coefficients in the asymptotic expansion: 

A4 



128 ' 
9a 2 - 48a 2 A 3 

M ~ 128 ' 

. 512^^ - 384^^3 - 9a 6 ^ . 
A ~ 1024 • (M) 

Next, we perform an asymptotic analysis beyond all orders to determine whether the 
solutions whose asymptotic behavior is given in (81) are stable. Let 

D(x)= yi (x)-y 2 (x), (85) 

where y\ and y 2 belong to the nth bunch of solutions. Then, D(x) satisfies the differential 
equation 

D'(x) = S[ndl - S[nd2 = sin ^ 1 ~ 6 ^ (86) 
cos 61 cos 6 2 cos 9\ cos 6 2 

If we now assume that D(x) is small, we obtain the asymptotic approximation 

X - 2 ~ {-Axy + 2ax)D{x) (x ->■ +00), (87) 

and since 

( ^ +00) - (88) 

we obtain the asymptotic approximation 

1 4 a 2 



cos 2 8 3 Ax 



2 



(x ->■ +00). (89) 
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Thus, D(x) satisfies the approximate linear differential equation 

D'{x) ~ D{x) + (x^ +00), (90) 

whose solution contains the arbitrary multiplicative constant K: 

D(X) ~ Ke l^m+a^/3 (a; ^ +00 ). (91) 

Note that as x — > 00, the difference function D(x) grows exponentially with increasing x, 
which contradicts the assumption that D(x) is small as x — > +00. We conclude that in 
the good wedge all solutions are unstable and that only a discrete set of isolated separatrix 
paths can continue deep into the good wedge without turning around and leaving the wedge. 
(For the case of the quantum harmonic oscillator the analogous unstable separatrix paths 
are shown in Fig. [7j) 

Numerical analysis (for the case a = 1) confirms that there exists just one unstable 
separatrix that runs directly from the left good Stokes' wedge to the right good Stokes' 



wedge. This path is displayed in Fig. 13 The path shown in this figure is the exact analog 
of the path in Fig. [7] that runs along the real axis from —00 to 00 for the case of the quantum 
harmonic oscillator. 

2. Asymptotic analysis of solutions in the bad Stokes' wedge above the positive-real axis. 



Next, we construct solutions to the differential equation (77) that approach the center of 
the bad wedge as x — > 00. Since the center of the wedge lies at an angle of ir/6 above the 
real axis, such solutions have the asymptotic form 

x „ Bo B-i Ba B^ Bp. , . , 

y-TM+^ + ^ + ^+^+^ + J--- (*->+«>)■ ( 92 ) 

\f rj 1K1 iXJ 1K1 iXj 1K1 

We determine B\ and B2 so that 9 approaches a constant as x — > +00. We find that if 

B 1 = \a and B 2 = fa 2 ^, (93) 

then 9 -> -4B 3 /y/3 as x 00. Substituting into the differential equation we get 1/V3 = 
tan(-4B 3 /V3), which gives 



73 / I 

— \ 



B 3 = ^(n--)7r. (9L) 



The higher-order coefficients are 

3V3a 4 

Bi ~ 
B 5 



128 

9a 2 - 48a 2 B? 



128 

-512^35! + 384X/353 + 9^^ 
B " ~ 1024 • (%) 

We perform an asymptotic analysis beyond all orders to determine if such solutions are 
stable. To do so, we define D(x) = y\{x) — y2(x), where y\ and 7/2 belong to the nth bunch 
of solutions, and we find that 
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FIG. 13: Separatrix paths for the ground state of the quasi-exactly-solvable VT anharmonic oscil- 



lator whose Hamiltonian is given in (24). (We have chosen the value of the parameter a to be 1.) 



The asymptotic behavior of these paths in the right good Stokes' wedge is given in (81). Exactly 
one separatrix path (heavy continuous line), which corresponds to the choice n = in (83), runs 
from the left good Stokes' wedge directly to the right good Stokes' wedge without veering into a 
bad Stokes' wedge. This path, which is the analog of the path along the real axis for the quantum 
harmonic oscillator, crosses the imaginary axis at y = —0.176 651 795 619 462 368. Other separatrix 
paths (heavy dashed lines) for the cases n = ±1 and n = ±2 are shown; these paths run from good 
Stokes' wedges to bad Stokes' wedges. The solutions to the differential equation are horizontal on 
the lightly dotted lines and vertical on the heavily dotted lines. 

where K is an arbitrary constant. Since D(x) vanishes exponentially for large x, the solutions 
in this bad wedge are all stable; that is, they bunch together as x — > +oo. The result here 



is analogous to that in (56) for the harmonic oscillator. 

Because these solutions are in a bad wedge, they blow up like the exponential of a cubic 
as they penetrate deeper into the wedge. Thus, the integral 

poo J 

Re [dzp(z)} = / ^^-2.73-^ w (97) 



diverges. Thus, the question is, Is it possible to have a contour that enters and leaves this 
wedge in such a way that this integral exists? That is, if y\(x) and y 2 (x) are two solutions 
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that enter that bad Stokes' wedge, does the integral 

poo poo / 2x 2 yi -2yf/3+ay 2 2x 2 y 2 -2y^/3+ay 2 

1= Re [dzp(z)] = / dxe-* ) (98) 

Jx=L Jx=L \ COS#i COS# 2 

converge? The answer to this question is no. 

The problem here is that the cubic polynomial in the exponent, 

Zx 2 y - \y 3 + ay 2 ~ ^^-x 3 + ax 2 + 0(x) (x ->■ oo), (99) 

causes the integrand of the integral to blow up like the exponential of a cubic. Thus, we 
cannot use the argument used earlier for the harmonic oscillator, where we were able to 
expand the exponentials and to calculate the difference in terms of hyperasymptotics. Thus, 
even if the integration contour enters and leaves the wedge, the resulting integral will not 
converge. Hence, it is not possible to leave this bad wedge as we did for the bad wedge of 
the harmonic oscillator; a contour that enters this bad wedge is trapped forever. 

3. Asymptotics of solutions in the bad Stokes' wedge centered about the negative-imaginary axis. 
This time we want to study the differential equation 

dx = cose 

dy sin 



and we want to investigate the integral in (78) for large negative y. This integral has the 
form 

p 2x 2 y— ax 2 

d3/e -V/3W£ (Id) 

sin# v ; 

A solution that approaches the center of the wedge has the form 

Cn n n n n 

x {y) — 2 + ^ + ~ + ~ + ^ + ~ + ••■ • ( 102 ) 
v y y v y y 

If we take 

C 2 = ad and C 3 = a 2 d, (103) 

then for large negative y we get 9 ~ —2C\ + O (y~ 3 ) and this balances dx/dy = O (y~ 3 ) if 
cos(— 2Ci) = 0, that is, if 

Ci = -|(n+|)7r. (104) 

We then get 

C 4 = ( G 3 - l)Ci, 
C 5 = i(2a 4 -5a)C 1; 

C 6 = i(2a 5 -9a 2 )d. (105) 
Thus, as y — > — oo, the asymptotic behavior of x(y) is given by 

, v „ (\ a a 2 a 3 -l 2a 4 - 5a 2a 5 - 9a 2 \ 
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Now we test for stability: Let 

D(y) = x x (y) - x 2 (y) 
be the difference of two solutions in the nth bunch. Then, 

cos#i cos 9 2 sin(6 l 2 — 6\) 



D'{y) 



sin 9\ sin 9 2 sin Q\ sin 9 2 



But for large negative y, 



9 2 -9 1 ~2D(y)(y 



2 2 

x 



ay) (y ->■ -00) 



The denominator is just [sin (n + |j 7r] 2 = 1, so (108) becomes 

D'{y)~2D{y){y 2 -x 2 -ay) {y -+ -00), 



and since x ^ C\y 2 , we get 



£>(y) 



2y 3 /3-ay 2 



2C\ 

3y 3 



(y ->■ -00) 



(107) 



108) 



(109) 



(110) 



(111) 



Thus, solutions that enter this bad wedge are stable as y — > —00. 

The crucial question is, Can the integration contour enter and leave this wedge? That is, 
can the contour go down the negative imaginary axis and then back up again. We use the 
integral in ( 78 ) to answer this question. We examine the integral 

^2x\y—ax\ ^lx\y—ax\ 



dy e 



dy e 



-2y A /3+aj/ 2 



sin Q\ sin 9 2 

. 2y 3/ 3+ay 2sm9 2 e 2x2y ~ ax2 — sm9 1 e 2x 2 y ~ ax 



_ Y sin 9 1 sin 9 2 

We approximate the denominator by observing that 9 ~ — 2C\ + O (y~ 3 ), so 



(112) 



sin^ ~ [(n + |) 7T + O (y- 3 )] 

~ (-ir C os[o(r 3 )] 

~ 1 + (r 6 ) (y^-00). 



(113) 



Thus, sin 9 ~ 1 + O (y ) as y — > —00, and we can replace the denominator in (112) by 1. 
To leading order the numerator in (112) is given approximately by 



sin # 2 - sintfi = 2 cos [\ {9x + 9 2 )] sin [\ (9 2 - 9 X 
~ 2 cos 9 sin [\{9 2 -9 1 )\ 

2A 



2 cos 



(n + §) 



7T 



sin [D(y) (y 2 



x 



ay)] 



sm 



2A 



■ ■ •) D(y) {y 2 



x 



ay) 



-4A(-l) n D(y)/y (y -+ -00). 



(114) 
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This gives a logarithmically divergent integral of the form J dy/y, just as we found in the 
case of the harmonic oscillator. 

We now follow the procedure that we used to analyze the harmonic oscillator. We examine 
the higher-order asymptotic behavior of the numerator and in expanding the exponentials, 
we include the first term beyond 1: 

sin 6 2 {2x\y — axf) — sin^i \2x\y — axl) . (115) 

We then add and subtract 



[sin 6*2 \2x\y — ax\ — 2x\y + aa^)] — [(sin#i — sin# 2 ) i^x\y — ax 2 ,)] • (H6) 

The second term in square brackets is negligible as y — > — oo. However, the first term is 
approximately 

sme 2 [(xl-xj)(2y-a)], (117) 
which leads to the asymptotic approximation 

OA 

{-l) n D{y) — {2y-2a)~4A(-l) n D(y)/y (y -oo). (118) 

yj. 

This exactly cancels the logarithmically divergent integral above and gives a convergent 
integral of the form J dy/y 2 . 



To summarize, we have shown that except for the unique path labeled n = in Fig. 13 



a path that solves the differential equation (77) [or equivalently, the differential equation 



(100)] in one of the good Stokes' wedges must leave the good Stokes' wedge and continue 



into one of the bad Stokes' wedges. As a result, the probability integral (78) [or equivalently, 



(79)] along such a path diverges exponentially. However, if the integral is taken along two 
paths, one that enters the bad Stokes' wedge centered about the negative imaginary axis 
and a second that leaves this Stokes' wedge, the integral along the combined path is finite. 



In Fig. 14 several complete paths that connect the left good Stokes' wedge to the right 



good Stokes' wedge are shown. There is a single path labeled (a) (this path also appears in 



Fig. 13 and is labeled n = 0) which runs directly from one good Stokes' wedge to the other. 
However, all other paths exhibit an intricate structure and repeatedly run off to and return 
from infinity in the bad Stokes' wedge that is centered about the negative-imaginary axis. 

B. Probability contours associated with excited states 

The eigenfunctions associated with higher energies have nodes in the complex plane. 
When J = 2, there is one node. The eigenfunctions have the general form 

^ n ( X ) = e -« 3 /3— 2 /2 (a; + c y (n9) 

The parameter c is determined by requiring that ip n (x) satisfy the time-independent 
Schrodinger equation 

d 2 \ 

x 4 + 2iax 3 + a 2 x 2 - Aix - E n %j) n [x) = 0. (120) 



dx 
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-3-2-10 1 2 3 



FIG. 14: Paths (heavy dashed lines) that connect the left good Stokes' wedge to the right 
good Stokes' wedge for the ground state of the quasi-exactly-solvable VT anharmonic oscilla- 



tor whose Hamiltonian is given in (24). There is exactly one such separatrix path (a) that 
goes directly from the left good wedge to the right good wedge, crossing the imaginary axis at 
y = —0.176 651795 619462 368. A path (b) that crosses the imaginary axis at a slightly higher 
point than the (a) path cannot reach oo in the good Stokes' wedge. Because they are unstable, 
such paths turn around, enter the upper bad Stokes' wedges, and can never re-emerge from these 
wedges. A separatrix path (c) is shown leaving the left good Stokes' wedge. This path enters 
the lower bad Stokes' wedge to the left of the imaginary axis, re-emerges along paths (d) or (h), 
and reenters the bad Stokes' wedge to the right of the imaginary axis. It then continues into the 
right good Stokes' wedge along the separatrix (e). Another separatrix path (f) leaves the left good 
Stokes' wedge and follows a more complicated course: After entering the lower bad Stokes' wedge, 
it leaves and returns along (g), leaves and reenters again along (d) or (h), leaves and reenters along 
(j), and finally enters the right good Stokes' wedge along the separatrix (f). Solution paths are 
horizontal on the lightly dotted lines and vertical on the heavily dotted lines. 

From this Schrodinger equation we obtain two equations, 3a — 2ic — E n = and (a — E n ) c = 
0, from which we conclude that the two energy levels and corresponding eigenf unctions are 

E = a ip {x) = e- ix3/3 ~ ax2/2 (x - ia), 

E x = 3a ^(x) = e~ ix3/3 - ax ^ 2 x. ( 121 ) 
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FIG. 15: Five paths (heavy solid lines) that emerge from the left good Stokes' wedge and attempt 
to reach the right good Stokes' wedge for the first excited state of the quasi-exactly-solvable VT 



anharmonic oscillator whose Hamiltonian is given in (24). The upper two paths veer upward into 



the bad Stokes' wedge and die there. The other three paths enter and reenter the lower bad Stokes' 
wedge and eventually pass through the node at the origin. Then, these paths repeat this process 
in the right-half plane and eventually succeed in reaching the right good Stokes' wedge. Note that 
there are six paths that enter the node at the origin; the two horizontal paths and the upper two 
paths eventually wind up in the upper bad Stokes' wedges. The lower two paths that leave the node 
enter the lower bad Stokes' wedge; these paths become part of the complicated route connecting 
the left to the right good Stokes' wedges. The solutions to the differential equation are horizontal 
on the lightly dotted lines and vertical on the heavily dotted lines. 



The probability contours associated with ipi(x) for the case a = 1 are shown in Fig. 15 This 



figure is quite similar in structure to Fig. 10 for the case of the quantum harmonic oscillator 



When J = 3, there are two nodes, and the associated probability contours associated 



with the highest-energy eigenfunction are shown in Fig. 16 The contours in this figure are 



qualitatively similar to those in Fig. 11 for the harmonic oscillator. 

As we move to higher energy and there are more nodes, the distribution of eigenpaths 
begins to resemble the canopy of the classical probability distribution in the complex plane. 
For the case J = 2 the classical canopy is shown in Fig. [17) 
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FIG. 16: Four paths (heavy solid lines) that emerge from the left good Stokes' wedge and attempt 
to reach the right good Stokes' wedge for the second excited state of the quasi-exactly-solvable VT 



anharmonic oscillator whose Hamiltonian is given in (24). The upper two paths veer upward and 



eventually curve around into the bad Stokes' wedge and die there. The other two paths enter and 
reenter the lower bad Stokes' wedge. After passing through both nodes these paths finally succeed 
in reaching the right good Stokes' wedge along separatrix paths. The solutions to the differential 
equation are horizontal on the lightly dotted lines and vertical on the heavily dotted lines. 

VII. FINAL REMARKS AND FUTURE RESEARCH 



In this paper we have shown that it is possible to extend the conventional probabilistic 
description of quantum mechanics into the complex domain. We have done so by construct- 
ing eigenpaths in the complex plane on which there is a real and positive probability density. 
When this probability density is integrated along an eigenpath, the total probability is found 
to be finite and normalizable to unity. 

There are many generalizations of this work that need to be investigated and much 
further analysis that needs to be done. To begin with, it is important to understand the 
time dependence of the complex probability contours. In this paper we have restricted 
our attention to the eigenpaths associated with eigenfunctions of the Hamiltonian. Such 
paths are time independent. However, for wave functions that are not eigenfunctions of 
the Hamiltonian, and even for simple finite linear combinations of eigenfunctions, there is a 
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FIG. 17: Analog of Fig. [4] for the quasi-exactly-solvable 'PT-symmetric anharmonic oscillator. The 
classical probability distribution in the complex plane for the case J = 2, a = 1, E = 3 is shown. 



complex probability current and the probability density flows in the complex plane. We have 
considered here only one elementary situation in which the complex probability contour is 
time dependent, and this was for the case of a complex random walk. We believe that a 
detailed study should be made of the high quantum-number limit; time- dependent contours 
C for non-eigenstates (such as gaussian wave packets) should be examined; the complex 
correspondence principle for coherent states should also be developed |50j . 

A second interesting topic for investigation is a detailed comparison of the classical pup 
tent in Fig. [4] and the analogous quantum picture. The pup tent in Fig. [4] was constructed 
by making the assumption that all elliptical paths in Fig. [3] were equally likely. But of course 
this is not quite valid. It is far less likely for a classical particle to be on a large ellipse than 
on a small ellipse close to the conventional oscillatory trajectory (the degenerate ellipse) on 
the real axis. A measure of the relative likelihood of being on any given classical ellipse is 
provided by the standard quantum probability density on the real axis as given in Figs. [T] 
and [2} Thus, we believe that the probability of being on a large ellipse is exponentially 
smaller than being on a small ellipse. With this improvement, the probability of a classical 
particle being somewhere in the complex- z plane can now be normalized to unity. (The 
volume under the pup tent in Fig. [4] is infinite.) 

With these changes in the distribution of classical probability, we can now begin to 
analyze the global distribution of quantum probability. The improved classical pup tent 
can now serve as a guide for estimating the relative probability of being on the various 
eigenpaths shown in Figs. [9 10 and 11 for the harmonic oscillator, and Figs. 14, 15, and 16 



for the quasi-exactly solvable anharmonic oscillator. We expect that the quantum pup tent 
describing the global distribution of quantum probability in the complex plane will have 
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ripples like the oscillations illustrated in Figs. [T] and [2] of the quantum probability density 
on the real axis. 

Finally, the most important feature of the quantum probability distribution in the com- 
plex plane - what we refer to above as the quantum pup tent - is that the density of 
probability along a complex contour is locally positive. However, the asymptotic analysis in 



Sees. [TV] and |VJJ of integration contours entering and leaving the bad Stokes' wedges leads 
us to conclude that for the total integral along an eigenpath to be finite some of the con- 
tribution to the probability integral must be negative. Our interpretation of this effect is 
not that the probability density is negative (the integrand is certainly positive), but rather 
that the contour goes in a negative direction and thus contributes negatively. (A trivial 
example of such behavior is given by the integral L dx. The area under the line y = 1 
is positive, but the integral is negative because it is taken in the negative direction.) This 
effect is interesting, and we believe that it deserves further examination. The fact that 
the total probability is unity but that individual contributions to the total probability are 
both positive and negative is strongly reminiscent of the results found in Ref. |5T] for the 
Lehmann weight functions for Green's functions of "PT-symmetric field theories. We be- 
lieve that the C operator, which is needed to understand the negative contributions to the 
Lehmann weight function, will ultimately play a significant role in the future analysis of the 
complex generalization of quantum probability. 
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